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We study the consequences of the f(R/0) gravity models for the Solar system and the large scale 
structure of the universe. The spherically symmetric solutions can be used to obtain bounds on the 
constant and the linear parts of the correction terms. The evolution of cosmological matter structures 
is shown to be governed by an effectively time dependent Newton's constant. We also analyze the 
propagation of the perturbation modes. Tensor and vector modes are only slightly modified, but 
two new scalar degrees of freedom are present. Their causality and stability is demonstrated, and 
their formal ghost conditions are related to a singularity of the cosmological background. In general, 
the Newtonian limit of these models has no apparent conflicts with observations but can provide 
useful constraints. 
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I. INTRODUCTION 

One covariant way of generalizing gravity is to consider the effective Newton's constant as a function of <f> = -R/D, 
the inverse d'Alembertian acting on the scalar curvature. This kind of nonlinear dependence can be parameterized as 
the following 

o : 

jag: j \ / 

General relativity is the limit / = 0. This form of the action can be motivated 1 by the fact that quantum effects 
£»\j , generically introduce corrections involving different operators and curvature invariants [l|, |2| • The exact form of the 
^ ' corrections is of course beyond control, but Eq.([T]) provides a simple parameterization, which can be constrained by 
observations: this is the phenomenological approach we take in the present study. 

Previously it has been noted that the most appealing (linear) form, f((f>) ~ <fi, could cure the unboundedness 
of the Euclidean gravity action [3|. More recently it was suggested Q that these types of effective gravity actions 
' . ■ could incorporate the observed cosmic acceleration with reduced fine tuning. The problem of the magnitude of the 
cosmological constant might be evaded since </> = R/d is a dimensionless combination, and thus no new scales (in 
addition to the Planck scale) need to be introduced. The coincidence problem could also get a simple resolution, 
since R and thus <f> vanish during radiation domination, so the acceleratory effect can be interpreted as a delayed 
consequence of the onset of the matter dominated era. 

Detailed analysis of the background dynamics shows that this could work, even for the simplest forms of / [f|. 
Unresolved issues exist though, like an occurrence of a sudden future singularity [a,ll|- In addition, in the simplest 
models featuring the late acceleration, one finds large effects taking place also during inflation. On the other hand, this 
could be used as beneficial effect to end inflation and reheat to universe, or even to drive inflation without the usual 
scalar fields as attempted in Ref. [a]. Realizing this possibility would require to find the classes of actions unifying 
consistently very early and the late universe. The fact that the action ([TJ) is capable of reproducing an arbitrary 
expansion history [5], hints that this is possible. However, in this work omit the background problematics, since our 
aim is to focus on inhomogeneous evolution of these models, particularly at the important Newtonian limit. 

This limit within the f(R) class of models has been intensively studied in connection with the dark matter and dark 
energy problems, see [3, [l(j, El \12, G3 for reviews. Even if viable examples would be found in that class of models, 
there seems to be no reason to prefer those more complicated models to the simple cosmogical constant, since they 
have not (at least this far) promised any help with the cosmological fine tuning problems. In this regard the f(R/D) 
scenarios make progress, as they, like mentioned above, could provide both a mechanism to connect the onset of the 
acceleration with the past matter dominated era, and remove the need for unnatural constants in the theory. Though 
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1 On the other hand, such models might appear in (stringy) compactification scenarios. This will become clearer when we introduce an 
equivalent description in terms of kinetically coupled scalar fields (non linear sigma model). 



also a f(R) action can be recasted into a scalar-tensor theory, the structure of the action (p}, which we'll soon clarify, 
is different. 

There has also been recently much interest in nonlocal scalar field cosmologies [1J|. These are motivated by the 
p-adic tachyon and open string field theory actions, which include operators with infinite number of derivatives acting 
on scalar fields. In particular, the exponent of the d'Alembertian, e D , appears. These kind of theories are equivalent 
to theories with several (possibly an infinite number of) coupled scalar fields, as soon becomes explicit in a different 
context. Nonperturbative methods can then be applied to treat localizable systems via higher-dimensional formulation 
HE EE Ell Ell • Friedmann- Robertson- Walker (FRW) solutions have been studied in string field theo ry Il9l . [20l |2ll |22| , 
and in related string-inspired forms of nonperturbative gravity with infinitely many derivatives [23l. l24j. 

Finally, let us remark that in the so called degravitation models the effective Newton's constant is promoted to 
a nonlocal operator, whose filtering effect could explain why gravity is insensitive to a cosmological constant [25| . 
leaving the possibility of the present cosmological term appearing as an afterglow of inflation [26( . These interesting 
models are burdened with problems associated to massive gravity [271 ] . including the loss of covariance, as the change 
of G is done at the level of field equations and not of the action. Though they are thus fundamentally different 
types of gravity modifications than we consider in the present paper, one can observe some analogous workings at a 
phenomenological level. In fact, similar actions could be constructed approximating the degravitation models 28]. 

The present paper is devoted to study of inhomogeneous solutions of the f(4>) case, and in particular in the weak 
field limit where the inhomogeneities appear as perturbations about more symmetric solutions. We study all types of 
cosmological perturbations for the first time, and look at the Solar system limit in more detail than previously done. 
In section (|TT| we present the model in generalized framework. In section ()III|) we write down the equations for each 
type of linear perturbations. The post-Newtonian limit, relevant for the Solar system physics and for the evolution of 
cosmological structure, is determined in section (|IV[) . The propagation of the scalar modes is more complicated than 
that of the vector and tensor modes, and we discuss that in the section (|Vj) . The conclusion is made in section (|VII) . 

II. THE MODEL AND ITS GENERALIZATIONS 

This section prepares for the study of perturbations. We present an equivalent description of action (fl"j) as a 
biscalar-tensor theory, that will be employed throughout the paper. We do this in a wider framework, looking first at 
more general operators than the d'Alembertian and then at more general functions / which may depend on several 
variables. Finally we write down the cosmological background system. 

A. Biscalar-tensor representation 

We consider the following class of nonlocal gravity actions, 

S=^Jd 4 x^dR(l + f(A- 1 R)). (2) 

where A is some operator of the form A = V^Q^V,,. A simple example is the d'Alembertian, □ = V^V^. A more 
complicated example is 

A P = □ + V^aR*" - pRg^Vv, (3) 

When a — 1 and (3 — 1/3, one has the Paneitz operator which is known to arise in conformal anomalies [291 ]. Einstein 
tensor -type model corresponds to (3 = a/2, and the box model is recovered when a = (3 = 0. Introducing the 
Lagrange multiplier £ to rename the inverse of A acting on R as 0, we rewrite the action in the form 

S=^Jd i x^—g[{l + f(4>))R + i{R-^)]. (4) 

Thus one may consider the theory as a multiscalar-tensor theory [301] . It is convenient to introduce a new field, 
;/> = .f{<P) — £, which is then the only one with a non-minimal gravity coupling. After a partial integration, the action 
assumes the form 

S= 2^J dAx ^~ 9 K 1 + ® R + Q ^ ( V ^ ~ /'(^) V m^) V„0 • (5) 



When this is coupled to the matter action, we get the field equations 
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where T™ is the matter energy momentum tensor. The effective energy momentum tensors for the two fields and for 
their coupling read as 
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Here S is a derivative operator, which includes derivatives up to m'th order, if Q^ v involves m'th order derivatives 

of the metric. More explicitly we can write 
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In the simplest case Q af} = g a/3 , we have T,^ af} = 5*6%. Only the sum of T* v + T$ v + T* is separately conserved. 
This implies the generalized Bianchi identity and usual conservation of the matter energy momentum tensor [311 ] . 
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B. More general models 

We consider a generalization of the model with an arbitrary number of derivatives. Let us then write 

s = ~h J dix ^9 R i 1 + f( A ~ mR > ■ ■ ■ . A ~ lR > AR > ■ ■ • > AnR )) • ( n ) 

This describes a large amount of possible models. The nature of these models is though not very transparent in 
this form. To better understand the field content and couplings of the model, we rewrite it in terms of scalar fields. 
Proceeding as previously, we arrive at the action 



S = i / d4x ^ K 1 + i>)R + Q^{V^ X i)(VvR) + V + K] 



(12) 



where the potential term is 
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and the kinetic term has the form 
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where / is a function of the n + m scalar fields, / = f(cj> m , ■ • • ■, 4>ii fit ■ • ■ > <Pn)- Nonlocality of the theory means that 
either m or n goes to infinity, and in a genuinely nonlocal one cannot reduce the number of derivatives (i.e. recast it 
in a local form); it is however common terminology to call e.g. the f(R/n) nonlocal. It is now clear that each higher 
power of derivative and inverse derivative in the action (|lll) adds two scalar degrees of freedom into the theory. At 
face value it seems that n + m + 1 of these are non-minimally coupled to gravity, however, we have shown that any 
nonzero n or any nonzero to both in essence add just one non-minimally coupled scalar into the play. In the Einstein 
frame the coupling resulting from nonzero to 7^ is transformed into non-minimal matter coupling. 

In passing, we remark that changing AR to MZ, where 1Z is some scalar quantity, would change only the ipR terms 
into iplZ terms in the actions ([5]) and (fl"2"|) . Any scalar 1Z which vanishes in de Sitter spaces might be useful, if one 
wanted to protect the de Sitter type solutions. 



It would be interesting to study in which cases the solutions may assume form of converging series. One of the 
interesting cases is the form / = logAi?, which represents the leading nonlocal loop correction to gravity as an 
effective quantum field theory (for previous approximative solutions, see [32, |33j). In the rest of the present paper, 
we will however specialize to the m = 1, n = case 2 . 

C. Cosmological equations 

Consider the inverse box model / = f(R/D). In a flat FRW background, 

ds 2 = -a 2 ( T )(dT 2 +dx 2 ), (15) 

we get the Friedmann equations 

3H 2 (1 + ^) = ^ + I (/'(^ _ ^) - 3#V> (16) 

- (2H + H 2 ){1 + i/>) = ^ + i (/'(0)0 2 - U) +i> + Hi, (17) 

Thus the time derivatives of the scalar fields contribute extra energy sources, and the field ip also modulates effective 
Newton's constant. Note that here dot means a derivative with respect ot the conformal time r, and H = (da/dr)/a 
is the conformal Hubble parameter. The equations of motion for the fields are then 

4> + 2HJ) = -6(H + H 2 ), (18) 

i> + 2Hijj = l"{m 2 - l2f'{cf>){H + H 2 ). (19) 

It will be useful to define the energy density fraction of matter as usual, 

n,„ = ^. (20) 

The background expansion of these models has been considered elsewhere [9, M, |34|. In particular, a power-law 
model / = fn4> n can feature modifications of the expansion both at early and at late late times, depending on the 
sign of (— l) n /„. Accelerating and super-accelerating scenarios exist both in vacuum and in presence of matter. An 
exponential model / = / e e A ^ has scaling solutions, thus the modification could be present at all times Q. In the 
following we will not specify any particular evolution but consider the general weak-field limit phenomenology of these 
models. 

III. PERTURBATIONS 

In this section we show the equations governing the evolution of cosmological perturbations. The vector and tensor 
perturbations are simple, but the analysis of the scalar perturbations, sourced both by the scalar fields and by matter, 
is continued in later sections. For cosmological perturbations in wide range of other modified gravity models, see [35j |. 

Introducing small perturbations about the background, one may write the metric as 

ds 2 = cl(t) 2 [-dr 2 (l + 2*) + drdx ■ W + dx t dx l (l - 2$ + H zj )] . (21) 

This is completely general parameterization of the inhomogeneities and anisotropics . Those can be decomposed into 
scalar, vector and tensor parts according the transformation properties under spatial rotations [3y]. At linear order, 
the different types of modes decouple. With the scalar perturbations, it is conventional to work in the Newtonian 
gauge [371 ] . That is specified by the metric potentials $ and 'J. Note that we use the uppercase symbols for the metric 



2 The generalization to the higher-derivative models is not completely straightforward. The potential l|13|l has to be taken into account 
whenever m > 1 or m > 0. If n remains zero, however, one then needs only to generalize all the (f> -dependent kinetic terms. 



potentials (and for all metric perturbations for consistency), which have nothing to do with the <\> and ip which label 
our scalar fields. The vector type perturbations are characterized by the W, which is transverse and thus has two 
independent components. The transverse and traceless 2?^ finally describes the gravitational wave. By construction, 
it is gauge invariant and includes two independent polarizations, which are conventionally called H + and Hx- In 
total, we then have 6 degrees of freedom in the perturbations, matching with the D(D + l)/2 — D physical degrees 
of freedom of the metric in dimension D — 4. The metric Eq. ([2Tj) is thus completely general and fixes the gauge 
uniquely as well. 



A. Scalar equations 

For completeness we list the equations for scalar perturbations here. The ADM energy constraint (G[j component 
of the field equation) is 
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and the momentum constraint (Gf component) is 
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The shear propagation equation (G* — \8 % Xi\ component) reads 
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where 7r m is the matter anisotropic stress. The Raychaudhuri equation (G\ — G® component) is now given by 



3$ + 3 2? 



1 



^ 



2(1 + V) 



$ + 3 H- 



1> 



2(1 + ib) 



* 



6H + -L-^ + 2f^-2^) + - 



* 



2 

2- (S Pm + 38p m ) + (2/^ - ipj 8<t> + f'^Scj) + 3^ - (3H + fySij>- (■ 



2(1 + V) 
The Klein-Gordon equation for the two fields read 

8'<p + 2H84> - V 2 8(f> = cj) (3$ + <$>) + (SH + 24> + Htf>) * - a 2 SR, 
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8ip + 2H8i>- \7 2 Sip = ip (3® + &) + (SH + 2i> + HiP + 2f" <f\ # + 2/"0<50 + (f'"<p 2 - 2f"a 2 R\ 8<f)-2fa 2 8R, (27) 
where the curvature scalar is 



R = 6(2? + H 2 ) + 2 -3$ - 9$ + 2V 2 $ - 3* - (62? + 6i? 2 + V 2 )* 



The matter components then obey 
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We consider later the forms of this system both in the vacuum and in the presence of matter. 



(30) 



B. Vector and tensor perturbations 

The find out how the rotational perturbation evolves, we need only the held equation 



ew=A r (p + p)— -v m , (31) 
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and the conservation equation 
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Combining these gives 

w+{2H+ ih> )w =T^w 1Tm - (33) 

Thus the rotational perturbation can grow only in the case dlog(l + ip)/dloga < —2. We can safely neglect the 
vector perturbations in these models, as usual in cosmology. 
The evolution of the gravitational waves is given by 

H + ,x + (™ + ^ H + , x + k*H + , x = ^^fZx, (34) 

so the different polarizations evolve in completely same way. Thus the waves still propagate with the speed of light 
and are sourced only by anistropic stresses of fluids. The decay of the amplitude of a cosmological background of 
gravity waves is given by the expansion rate in a way which is only slightly modified by the evolution in tp. 

The propagation of scalar modes in these models turns out to be more involved, and will be considered in the 
separate section (|V|) . 

IV. NEWTONIAN LIMIT 

A. Solar system 

Consider first the Newtonian limit of the cosmological perturbations to study the effects within our Solar system. 
In the Solar system we may neglect the cosmological expansion and set a(r) = 1. Furthermore, one may assume the 
time derivatives are negligibly small with respect to space derivatives. We consider vacuum outside a mass and thus 
can neglect the source terms. The ADM energy constraint (j2"2")l then reduces to 

V 2 $ = ; 1 .. V 2 ^. (35) 

2(1+ V) 

Assuming spherical symmetry, the solution to this equation is $ = Sip + C\/r + CV Demanding that we recover 
the usual Schwarzschild solution in the limit of GR fixes the two integration constants as C\ — —2GM, C2 = 0, 
where M is the mass of the spherical object. This shows that the gradient of the held ip ac t s as an extra source in 
the post-Newtonian generalization of the Poisson equation. The shear constraint (|24[) shows that the gravitational 
potentials are not equal due to the effect of this held, 

1 + tp 

The Raychaudhuri equation ()25|1 is redundant with the above constraints at this limit. We then exploit the Klein- 
Gordon equation (|2T|) to obtain the constraint 

6i/> = -4/'(0) (* - 2$) , (37) 

where we have integrated again twice and killed the integration constants by requiring consistency at the GR limit. 
Now we can eliminate the gradient 5ip from the above equations. The solution for the gravitational potentials is then 

x5j = _2_GM/l ± iP-8f\ $ = _2GM/ l + ^-4/' \ 

l + j> \l + ip-6f'J ' ' l + ip\l + ip-6f'J' v ; 



where the background values ip and /'(</>) are constants. So, we may write the metric as 

ds 2 = - ( 1 - *™£) dr> + (l + *^ dr 2 + rW, (39) 



where the effective gravitational constant is 
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and the post-Newtonian parameter 7 is given as 



7 Vi + ^-8/7 

There are very tig ht constraints on the time variation of G* and the deviation of 7 from its value in general rel- 
ativity 7=1 [381 ] . As the current accuracy in determination reaches few parts in hundred thousand, the most 
stringent constrain ensuing the tracking of the Cassini spacecraft [39f, we can estimate on that our fields should 
satisfy \f'(<j))\, \ip\ ^ 10~ 4 . If we would naively set the field <j) to its cosmological value in the models where the /(</>) 
correction drives the present acceleration, those models would be immediately ruled out! However, since the exact 
Schwarzschild solution is Ricci-flat, R = 0, for a properly regularized operator one should have <fi = d^ 1 R = there. 
Constraints can still be imposed. For this purpose, expand the coupling f(<j>) as a power series, 

00 

/(fl = X>^- ( 42 ) 
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Negative powers are excluded since we want Ricci-flat spaces to exist. We find that 

f l + fo-8.f 1 \ G /1 + jWM ,.„. 

Therefore it is possible to constrain the constant and the linear parts of the coupling from Solar system experiments. 
However, the constant /o may always be absorbed into a redefinition of the k — > k/(1 + /o) in the action (fT]), and so 
one may set /o = 0. We then get the strict bound 

-5.8-1CT 6 < /1 <5.7-l(T 6 (44) 

Deviations from GR ensuing from the higher order terms in the expansion for /(</>) are of the second order in 
perturbation theory, ~ (GM/r) 2 . Since the gravitational potential $ = —GMq/tq ~ 10 -6 at the surface of the Sun, 
we expect the order of magnitude of the nonlocal corrections of the type /2</> 2 + fs4> 3 + ■ ■ ■ in these models to be within 
the experimental bounds based on the Schwarzschild geometry. 

B. Structure formation 

Next we consider cosmological scales and relax some of the simplifying assumptions. We then allow time evolution 
and do not impose spherical symmetry. We take into account pressureless matter source and thus focus on the 
formation of inhomogeneous structure in dust-dominated cosmology. We look at the first the subhorizon limit. Then 
the gradient terms are more important than the time derivatives of the inhomogeneities, y ~ H 2 y -c k 2 y. This 
approximation could be called the linear post-Newtonian limit, as in [401 ] . 

Differentiating Eg. (f!?9"]) and using Eq. (|3"0]) quickly gives us 

o m + H5 m = kH. (45) 

We need to solve the right hand side of this equation. The response of the gravitational potential to the matter source 
is now modulated by the gravity coupling, 

^ = 5oTvo(^ + ^). (46) 



The Raychaudhuri Eq. (125)) at this limit gives 

^ 2 ^_l_(^ m + fc ^). (47) 

It is through these two relations that the modified gravity effects could now show up in the integrated Sachs- Wolfe 
effect and in the weak lensing potential. These effects and their correlation with the cosmic microwave background 
anisotropics could thus be used to probe these gravity models. The latter is relevant at the large scales, where the 
use observations is limited by the cosmic variance. From here we see that the relative change of the potentials is of 
roughly of the order -0, which is typically negative and order of 0(0. 1) in the simplest models accelerating the universe 
nowadays [5j. This is in a viable but interesting range. 

Combining Eqs. ()46|47|) with Eq.((27|) gives us the gradient of ijj in terms of 5p m - Then Eq. (|45|) becomes 

6 + HS = AnG^ Pm S = ( Hilly ) TT^ PmS - (48) 

where the expression for the effective gravitational constant is the same as in Eq. (|4"0"|) . Now it is, unlike in the near- 
Schwarzschild approximation, in general time-dependent, if the background <p and ip are evolving. A crucial feature 
here is the absence of effective sound speed terms. These are known to appear in some coupled dark energy [4J| and in 
modified gravity [42], causing then instability in the dynamics of linear perturbations. This can be tightly constrained 
by observations of the galaxy distribution. However, our result (|48| suggests that the growth of structure is here 
scale- independent, and thus the constraints will now ensue solely from the overall normalization and the growth rate. 
Since the G^/G is of order one or less in the /(</>) dark energy models, we expect of at most of order one modifications 
to the growth rate at late times. This is within the observed limits at the present [43(, but could be a useful additional 
constraint on these models. 

Finally we remark that all of the results about the large-scale structure formation in this subsection rely crucially 
on the assumption that dark matter is exactly cold (a perfect and pressuless fluid). By relaxing this assumption one 
may drastically change the predictions of modified gravity models concerning the inhomogeneities, while keeping both 



the gravity sector and the background expansion fixed [44 1 



V. GHOSTS, CAUSALITY AND STABILITY 

To study the propagation of the scalar modes, it is useful to go to the Einstein frame. Perform the Weyl rescaling 

mm 

g^ = (l + ^)T^g^ = eT^ x g^. (49) 

For a moment we consider the D-dimensional case. The action then becomes [5| 
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Now %j} has also acquired a kinetic term. The curvature coupling has been removed in this frame, but the matter 
sector has now become nonminimal. However, this does not complicate our analysis if solely the vacuum is considered. 
Furthermore, it is possible to to diagonalize the kinetic Lagrangian to get rid of the derivative interaction. For that 
purpose, we must define the new field £ 
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The diagonalized action is then in vacuum 
d B x-s[-J) 
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In presence of (non-conformal) matter both fields would now be non-minimally coupled to the matter sector. One 
notes already that in the / — > limit we have one canonic and one phantom scalar field (except in D = 2 where the 
phantom field disappears) . Indeed we will find that this is the usual field content of these models also in the more 



general case. In the rest of this section, we omit the tildes from the Einstein frame quantities to ease notation, since 
every variable will be evaluated in the Einstein frame except when explicitly mentioned. 

When analyzing perturbations in multi-field systems it is convenient to perform a decomposition of the perturbations 
into adiabatic and entropic parts [47| . In our two-field case the entropy will then be presented by just one field. Except 
for possible differences in signs, our form ([52]) belongs to the non-linear sigma models studied in [48(. For a recent 
consideration on more general non- linear sigma models, see [49 ]. To proceed, set now D = 4 and call 



F(60 = eVft- 

The background adiabatic field a is implicitly given by 
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and the entropy field is orthogonal to this in the field space. The entropic and adiabatic perturbations can then be 
written as 






(55) 

(56) 

With the four previous definitions, the Klein-Gordon equations for the fields </> and £ and a lot a algebra, one obtains 
the simple results 



5a + 2HSa + k 2 Sa = & (a - 3Ha + an) 



(57) 
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6s = 0. 



(58) 



Thus we have found the two independently evolving scalar degrees of freedom of these models. Coupling of the 
evolution equations seems to be a almost generic property of generalized multi-field models [4!| ■ The fact that the 
adiabatic and entropy field evolutions now decouple is due to the masslessness of the fields (they do not have potentials, 
self or otherwise), and the diagonalizability of the kinetic Lagrangian in the Einstein frame. One may choose a gauge 
where the left hand side of Eq. (|57p vanishes, but in general the perturbations of the a- field are sourced by the metric 
perturbations. The evolution equation for the entropy mode, on the contrary, is homogeneous. This means that if 
the s-ficld is initially smooth, it will stay so. At small scales any initial perturbation will be washed away because of 
the sound speed equals unity, but at large scale there if a potential instability if \[F \^/ ' \f~F~ > 0. Most importantly, 
also the other field has the propagation speed c A equal to the light speed, just like it is for the canonical scalar fields, 
c A = 1. Therefore, in all directions of the field space the propagation speed is unity. Since the Weyl rescaling retains 
the light cones, the scalar perturbations propagate in both frames with the speed of light. Thus we can a posteriori 
justify performing the analysis in the Einstein in this case, though in general the propagation velocity of perturbations 
may not be frame-invariant. 

Thus, at the classical level there is no instability (since c A > 0) or causality violation (since c A < 1). At the 
quantum level, an inconsistency may occur if the kinetic terms of the fields have a wrong sign. Our derivation 
above assumed that F(<fi, £) > 0. If this is not case, the kinetic term of <j) would have the wrong sign, implying an 
appearance of a ghost at the quantum level. If furthermore _F(</>,£) < (£/0) 2 , the adiabatic field direction in the field 
space corresponds to a ghost degree of freedom. Translated into a condition for the original Jordan frame fields, the 
no-ghost requirement for the scalar modes reads simply 



1 



1 + ip 6 



>0^>6f{(f>) > 1 + V, 



(59) 



where in the second line we have used the tensor no-ghost requirement. In fact, the effective equation of state of the 
(Jordan frame) universe can be expressed as 



W eff 



3 Wm n m + rmh-kii 2 - 6) + rw^ 2 - (i - m^M 

3(l+V-6/'(«£)) 



(60) 
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Note that this is fully general expression and does not assume vacuum. Thus the future sudden singularity occurs 
just when the limit ([55]) would be crossed. Is then clear that kinetic term of the scalar mode does not change its 
sign. However, we can make the observation that exactly when the inequality ()59|1 becomes an equality i.e. the scalar 
kinetic term changes it's sign, the gravitational constant P0|) diverges. In fact this applies also for the tensor mode 
kinetic term, since G* ~ 1/(1 + ip)- The same phenomena has been observed to occur in Gauss-Bonnet cosmology. 
There exist peculiar instabilities, of stringy vacua [501 ] and of linear tensor perturbations [5lL |52|. In particular, a 
divergence of scalar modes has been associated with a crossing of a ghost limit [53j, |54I |. There, in contrast to our 
present case, the background can continue to evolve smoothly while the singularity occurs only at the linear level. 

VI. CONCLUSIONS 

We considered a class of cosmological models based on the generalized gravity action |T]). Though apparently 
nonlocal, this action may be recasted into a local biscalar-tensor theory. The theories represent a subclass of more 
general covariant theories, involving more, possibly infinite number of scalar fields. The simple parameterization of the 
corrections to Einstein gravity as the function f(R/U) has already several new theoretical features and potential to 
describe new phenomenology. In particular, one may construct dark energy cosmologies without introducing unnatural 
parameters and explain the coincidence of the similar matter and dark energy densities just today by the response of 
R/O to the onset of matter domination. 

In this paper we focused on the cosmological perturbations and weak field limit of these models. The scalar, vector 
and tensor perturbation equations were derived up to linear order in the flat FRW background. To simplify the 
equations, we fixed a gauge for each type of perturbation, but of course the equations in an arbitrary gauge may be 
obtained from our results using the standard transformation formulas. The scalar equations we wrote down in the 
Newtonian (or longitudinal) gauge. We considered the evolution of inhomogeneities in a matter dominated universe. 
We found that gradient effects are absent, but that the growth rate is modified according to a time evolving effective 
gravitational constant. Also, as rather generic predictions of modified gravity, the Poisson equation is modified and 
there is effective anisotropic stress. These effects might provide hints to detect these models from any observed 
deviation from the general relativistic predictions. In particular, the integrated Sachs-Wolfe effect probes the late 
time variation of the gravitional potentials [5a . l5q | , and the small scale effects of modified gravity might be probed by 
weak lensing experiments [571 158| | . We hope to quantify these constraints on specific models by a detailed numerical 
analysis in a future work. In addition, inflation and the generation of the primordial spectrum would be interesting 
to study with the action {1}. 

We considered here the Post-Newtonian limit for static spherically symmetric solutions, which is crucial for the 
local tests of general relativity. We found that there are two PPN type parameters possibly differing from their 
general relativistic values, with a given dependence of the form of /. This dependence allows to constrain tightly 
the constant and the linear parts of the function /. In particular, they are severely restricted by the bounds on the 
post-Newtonian 7 obtained from the tracking of the Cassini spacecraft, resulting in Eq. (144|) . For example, the purely 
linear model resulting in Q m = 0.3 today is ruled out by this constraint alone (though this model would not accelerate 
the universe enough anyway Q), as well as the simplest linear stabilator of the Euclidean action (though it already 
had tension with the Big Bang nucleosynthesis (3j). The constraints on the higher powers in / seem to be some orders 
of magnitude looser. 

The issues of causality and stability were also considered. The rotational perturbations were shown to decay as 
usually. The gravitational waves evolve with an extra friction term and a modulated response to the matter stress 
sources due to the coupling to (1 + ip). The evolution of the two extra scalar degrees of freedom is more delicate. 
We found that it is possible to separate these degrees of freedom into decoupled modes by considering the theory 
in the Einstein frame, diagonalizing the kinetic terms and then considering the perturbations Ss and da, which 
are fluctuations of implicitly defined combinations of the fields appearing in the original Jordan frame Lagrangian. 
However, since it is possible to write down decoupled, canonical actions for those two perturbations, we believe it is 
legitimate to regard them as the physically propagating degrees of freedom. This recipe would not probably work 
in a more general case though. Here our result is that the new scalar modes propagate perturbations with the 
light speed, and thus they are causal and stable. The formal ghost conditions are directly linked to divergence of 
classical perturbations (similar to the Gauss-Bonnet cosmology) and here also with the sudden future singularity of 
the cosmological background. Thus these conditions clearly indicate a pathological point of the solutions. We note 
that the singularity structure of the f(R) models has also raised recent concerns [71 [53. [60L l6l| . 

To conclude, the f(R/U) gravities seem have a viable Newtonian limit, although they may feature significant 
modifications of the cosmological expansion. 
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